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Abstract 

CNj , The inherently homogeneous stationary-state and time-dependent Schrodinger equations are often recast 

into inhomogeneous form in order to resolve their solution nonuniqueness. The inhomogeneous term 
can impose an initial condition or, for scattering, the preferred permitted asymptotic behavior. For 
bound states it provides sufficient focus to exclude all but one of the homogeneous version's solutions. 
■ Because of their unique solutions, such inhomogeneous versions of Schrodinger equations have long been 

the indispensable basis for a solution scheme of successive perturbational corrections which are anchored 
, by their inhomogeneous term. Here it is noted that every such perturbational solution scheme for an 

inhomogeneous linear vector equation spins off a nonperturbational continued-fraction scheme. Unlike 
its representation-independent antecedent, the spin-off scheme only works in representations where all 
components of the equation's inhomogeneous term are nonzero. But that requirement seems to confer 
theoretical physics robustness heretofore unknown: for quantum fields the order of the perturbation places 
i <~| a bound on unperturbed particle number, the spin-off scheme contrariwise has only basis elements of 

Cl^i unbounded unperturbed particle number. It furthermore is difficult to visualize such a continued-fraction 

^J_| ■ spin-off scheme generating infinities, since its successive iterations always go into denominators. 

Introduction 

Schrodinger equations, whether stationary-state or time-dependent, are homogeneous, and as such can suffer 
from solution nonuniqueness. In order to be able to effectively apply standard successive approximation 
schemes, such as perturbational ones, to Schrodinger equations, they are often first recast into a specialized 
inhomogeneous linear form which has a unique solution. Denoting its state vector as \ip), a simple generic 
presentation of a Schrodinger equation in inhomogeneous form is, 



|^) = ^o)-V|V). (la) 

> 

The state vector l^o), which comprises the inhomogeneous term of Eq. (la) usually shares a key feature with 
\ip), such as asymptotic behavior in a region of infinite extent or exact value at an initial time, or, alternately, 
is a passable approximation to \tp). Given that the Hamiltonian operator of the homogeneous Schrodinger 
equation which underlies Eq. (la) is H, we in addition assume the existence of an exactly diagonalizable 
Hamiltonian operator Hq that has \ipo) as one of its eigenstates , and for which the linear operator V of 

C*") | Eq. (la) is of no less than first order in V = (H — Hq). 

If we now take the viewpoint that the inhomogeneous term l^o) on the right-hand side of Eq. (la) is a 
sufficiently good approximation to on its left-hand side that the remainder term —V\ip) is only a small 
■ perturbation, then the perturbational successive approximations, 

= IV^o) - V|V (n) ), (lb) 

where, of course, \tp^) — \tpo), would be expected to converge reasonably rapidly. 

Additional insight into the Eq. (lb) perturbational iteration scheme can be obtained from taking the 
second term on the right-hand side of Eq. (la) to its left-hand side, which yields, 

|V)+V|V) = |Vo). (1c) 

Eq. (lc) has the formal operator solution, 

|V) = [l + V]- 1 |^o). (Id) 
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This presents the issue of how to practically evaluate the formal operator [1 + V] 1 . One way is to try to 
use its formal geometric series expansion, 

[1 + V]- 1 = 1 - V + (V) 2 + • • • + (-V)" + • • • , (lc) 

which is readily verified to yield the exactly the same result when inserted into Eq. (Id) as is obtained from 
successively applying the Eq. (lb) perturbational iteration scheme. Given the well-known convergence and 
divergence characteristics of the geometric series, it is apparent that the Eq. (lb) perturbational iteration 
scheme for the inhomogeneous Eq. (la) is unlikely to be greatly useful unless \ipo) * s substantially dominant 
over V | -00 ) ■ 

That naturally raises the question of whether there might be an alternate general iteration approach 
to the inhomogeneous Eq. (la) that could come to the rescue when its Eq. (lb) perturbational iterations 
converge too slowly or diverge. A Lippmann-Schwinger nonrelativistic potential-scattering variant of the 
inhomogeneous Eq. (la), considered for arbitrarily strong potentials in a recent article [1], has revealed 
that if a complete representation basis set can be found such that every component (pi\4>o) of the 

inhomogeneous-term state vector |?/>o} is nonzero, then a robust nonperturbational iteration scheme can be 
devised from the explicit presentation of Eq. (lc) in that representation, namely from the specific equalities, 

{p^} + ( Pi \V\^) = (pM- (2a) 



If we now make the assumption that (pi\ip) is, like (pi\ipo), nonzero, we can factor the left-hand side of 
Eq. (2a) into (pi\ifj) and the resulting nonzero cof actor. We then proceed to divide both sides of Eq. (2a) by 
that nonzero cofactor, which changes the appearance of Eq. (2a) to, 

{PiM = (^ivo)/ [i + ((p^r 1 (p»iviv)] , (2b) 

from which we straightforwardly devise the manifestly nonperturbational successive-approximation scheme, 

= <P*lV>o>/ [l + ((P^W 1 (Pi\V\^ n) )} , (2c) 

where, of course, (pi\ip^) = {pi\ipo}. If ( j Oi|^"- ) ) is nonzero, then because (pi\ipo) is nonzero, barring the 
pathological occurrence of a divergence in [pi\V\ip^) , {pi\^ n+1 ^) will in turn be nonzero. We also note that 
the iteration scheme of Eq. (2c) has the desirable nonperturbational character of a continued fraction. 

Thus the Eq. (la) generic specialized inhomogeneous form of the Schrddinger equation always spawns 
not only the perturbational successive-approximation scheme of Eq. (lb), but as well the nonperturbational 
continued-fraction successive-approximation scheme of Eq. (2c) that can regarded as its spin-off. 

We shall now survey some of the well-known circumstances where a Schrodinger equation with Hamil- 
tonian operator H for a state is combined both with a state \ipo) which has a crucial similarity to \ip) 
and as well with an exactly diagonalizable Hamiltonian operator H a that includes \tp a ) as one of its eigen- 
states to produce an inhomogeneous linear equation for \ip) which has the generic form given by Eq. (la), 

^ ~~ def ^' ^' 

where the linear operator V is of at least first order in V = (H — H ). We shall as well mention some of 
the complete representation basis sets {{pi\} for which certain of the surveyed inhomogeneous terms | Vo) 
have exclusively nonzero components (pi\tpo) — note that the eigenstate set of the Hamiltonian operator H 
is entirely unacceptable as such a complete representation basis set {{pi\} because all of its members aside 
from (tpo\ itself are orthogonal to |-0o) - 



The Lippmann-Schwinger equation for nonrelativistic potential scattering 

As in Ref. [1] we consider the coordinate-representation nonrelativistic Schrodinger equation for an eigenstate 
(r\ipE) of positive energy E, 

(-ft 2 V 2 /(2m) + V(r)) (r\^ E ) = £(r# S ), (3a) 
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where, 

lim V(r) = 0. (3b ) 

|r|— ¥00 y ' 

Thus for sufficiently large |r|, the Schrodinger equation of Eq. (3a) reduces to, 

(-h 2 Vl/(2m)) (t\1> e ) = E(r\i, E ), (3c) 

which is satisfied by any plane wave e tp ' r l n for which |p| = (2mE)i and by any linear superposition of 
these plane waves. Among those linear superpositions are all the angularly modulated ingoing and outgoing 
spherical waves that have wave number k = (2m£)J /%. 

Now a scattering experiment at energy E > is described by a particular solution of the Schrodinger 
Eq. (3a) which at sufficiently large |r| that Eq. (3a) is well-approximated by Eq. (3c) consists of only a single 
specified plane wave e lpT l h of momentum p plus only outgoing spherical waves [2]. We denote this scattering 
solution of the Schrodinger Eq. (3a) as (r| ?/;+). 

It turns out that an inhomogeneous modification of the Schrodinger Eq. (3a) describes (r|i/>+) uniquely, 
namely the following nonrelativistic Lippmann-Schwinger equation for potential scattering [2], 

<i#+) = e***'* - {t\{H q -E p - ie)- 1 ^), (3d) 

where H = f |p| 2 /(2m) = — ?i 2 V 2 /(2m) is the kinetic energy operator and E p = f p| 2 /(2m) is the kinetic 
energy c-number scalar that corresponds to the c-number momentum vector p. 
Taking (r|p) = e' p ' r / R , we can write Eq. (3d) in the form of Eq. (fa), i.e., 

\i>+) = \p)-(H -E p -ie)- 1 V\i>+), (3e) 

where |-0o) = |p) and V = (Hq — E p — ze) -1 ^. Furthermore, since Ho\p) = E p \p), we can recover the Schro- 
dinger Eq. (3a) from Eq. (3e) by multiplying the latter through by (H — E p ), followed by rearrangement of 
the resulting terms between the left-hand and right-hand sides. 

The negative imaginary infinitesimal — ie that appears in the Lippmann-Schwinger Eqs. (3d) and (3e) 
ensures that only outgoing spherical waves are present for sufficiently large |r|, in addition, of course, to the 
single specified plane wave e tp ' r / h of momentum p. 

The usual perturbational iteration of the Lippmann-Schwinger Eq. (3d) is, 

(r|4" +1)+ ) = e'P-'/ft - (r\(H -E p - ie)-'V\^ p n)+ ), (3f) 

which with {r\ip p ^ + ) = e tp ' r / h generates the familiar perturbational geometric Born series for nonrelativistic 
potential scattering [3]. 

For nonperturbational "continued-fraction" iteration of the Lippmann-Schwinger Eq. (3d) we happen 
to be in the extraordinarily fortunate situation that coordinate representation of the inhomogeneous term 
IV'o) — |p) is always nonzero because (r|p) = e lp ' r l h ^ 0. Therefore in the manner of Eq. (fc) and Eqs. (2a)- 
(2c), the Lippmann-Schwinger Eq. (3d) in coordinate representation gives rise to the nonperturbational 
"continued fraction" iteration scheme, 



f + (<r|4"> + >) 1 (r\(H -E p - ie)~'V\^ + ) 



(3g) 



where we of course have that (v^p ^) = e lp ' r / n . 



A bound-state inhomogeneous equation suitable for iteration solution 

The just-discussed Lippmann-Schwinger equation is an inhomogeneous Schrodinger-equation variant that is 
of the form of our generic Eq. (fa). With an arbitrary plane wave as its inhomogeneous term; the Lippmann- 
Schwinger equation enables completely prescribed iteration refinement of that approximating plane wave 
toward its unique associated "asymptotically outgoing-only spherical-wave" exact scattering wave function. 
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The key ingredient that is needed to fashion the inhomogeneous Lippmann-Schwinger equation from its 
underlying homogeneous Schrodinger equation is of course the "outgoing-only spherical-wave" free-particle 
propagator, which is constructed from the selfsame free-particle Hamiltonian for which that approximating 
plane wave is an eigenstate. 

In contrast, the standard Schrodinger approach to the perturbational refinement of an approximating 
bound state \ipj) fails to explicitly present the inhomogeneous Schrodinger-equation variant which links that 
approximating bound state to a unique associated exact bound-state solution of the underlying homo- 
geneous stationary-state Schrodinger equation H\tpj) = Ej\ipj). Due to the absence of the information which 
is inherent to that inhomogeneous linking equation, the standard Schrodinger bound-state perturbational 
approach can't offer minutely prescribed iteration instructions for its perturbational refinement procedure [4] . 
More importantly, if that inhomogeneous linking equation is not in hand, there is no obvious way to devise 
alternate nonperturbational methods for the iteration refinement of approximating bound states. 

Fortunately, we can follow the blueprint that is provided to us by the Lippmann-Schwinger equation to 
devise its missing bound-state analog. A major difference between the two, however, is that bound states 
certainly don't feature outgoing, ingoing or any other type of traveling wave. Therefore the kind of propagator 
we need for bound states is of the standing-wave type. Nor can a propagator suited to an approximating 
bound state possibly be constructed from the free-particle Hamiltonian operator as it properly is for the 
approximating plane waves of the Lippmann-Schwinger equation. The construction of the standing wave 
propagator we need must be from a Hamiltonian operator Ho which has that approximating bound state 
as one of its bound-state eigenstates, i.e., H must satisfy, 

H \^)=E^). (4a) 

Of course, since that approximating bound state |V>°) is indeed bound, we can normalize it to unity, i.e., 

(V>X> = 1. (4b) 

Finally, the construction of the required standing-wave propagator from the Hamiltonian operator H can't 
be accomplished in practice unless H n is exactly diagonalizable. If all the eigenstates and eigenvalues of H 
are actually available, its corresponding standing- wave propagator can be obtained in the schematic form, 



P/(H E°)} ^ lim (H Ef) Uh E°) 



(4c) 



where the letter P in the standing- wave propagator's definition denotes "principal value". The exactly 
diagonalizable Hamiltonian operator Hq is frequently called the "unperturbed" Hamiltonian operator, its 
particular bound eigenstate |^) the "unperturbed" approximating bound state, and that state's Hq eigen- 
value Ej the corresponding "unperturbed" approximating bound-state energy. By the same token, one could 



call the operator 



P/{H -®) 



P/(Ho-Ef)] 



the "unperturbed" standing-wave propagator. Note that 

has a special feature which has no analog for traveling-wave propagators such as the outgoing-wave propa- 
gator of the Lippmann-Schwinger equation, namely that, 



m p/(h -e") 



= 0. 



(4d) 



The final ingredient we need to obtain the bound-state analog of the inhomogeneous Lippmann-Schwinger 
equation is the homogeneous stationary-state Schrodinger equation for the exact bound state \tpj) that is 
uniquely linked to the "unperturbed" approximating bound state | ?/>"). Denoting the exact Hamiltonian 
operator for the physical system we are studying as H, this homogeneous stationary-state Schrodinger 
equation is of course, 

H\i> j )=E j \i> j ), 
which it can sometimes be more convenient to express as, 



H-E; 



\1>i)=0. 
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Since the physical system's exact Hamiltonian operator H won't in general be exactly diagonalizable (unlike 
the somewhat artificial "unperturbed" Hamiltonian operator Ho), we have no more a priori knowledge of the 
particular exact bound-state energy Ej than we do of its corresponding exact bound-state eigenstate \ipj) 
of the exact Hamiltonian operator H . In order to avoid having the exact bound-state energy Ej become a 
completely independent unknown for which we need to solve, we choose to write the physical system's exact 
homogeneous stationary-state Schrodinger equation for the bound state \tpj) in the more cumbersome form, 



H 



0. 



(4e) 



which entirely supplants Ej by \ipj) and H. 

Eqs. (4a)-(4e) supply all the ingredients needed to devise the bound-state analog of the Lippmann- 
Schwinger equation. Inspection of the Lippmann-Schwinger equation strongly suggests that we multiply 
the physical system's homogeneous Schrodinger Eq. (4e) from the left by the "unperturbed" standing-wave 



propagator 



P/(H — Ej) as the first step toward that analog. That operator multiplication produces the 



physical Schrodinger-equation variant, 

~P/{H -E°j) 



H 



0, 



(4f) 



which is not yet the inhomogeneous Lippmann-Schwinger equation analog that we seek. To complete our 



to which we shall then add 



task we now develop an identity involving \tpj), [p/(Ho — Ej) and — Ho + Ej 
the Eq. (4f) physical Schrodinger-equation variant. One might naively expect that, 

\4>j) + [P/(H - Ef)] [-Ho + E%] \i>j) = 0, 

but if we multiply the left-hand side of this proposed equality by {ipj | and take note of Eq. (4d), we see that 
the result is (i>j\ipj) rather than zero. From Eqs. (4c) and (4b) it is indeed seen that, 

IV;) + [P/(H - £$)] [-H + 25?] |^) - (4g) 
We now add the Eq. (4g) identity to the physical Schrodinger-equation variant given by Eq. (4f) to obtain, 



Wj)+ P/(H -Ef] 



H - Hi 



-0 + 23? - (<Vil%i>/<m->)] Wi) = W^Mi)- 



If we now define the "interaction Hamiltonian operator" V as, 



V=(H-H ), 



we can rewrite Eq. (4h) as, 



|^> + [P/(H Ej)] [V ^j\(V + H - E^j)/(^j))] \^j) = 



(4h) 



(4i) 



(4j) 



If we move the term involving the "unperturbed" standing- wave propagator 
side of Eq. (4j), we obtain, 



P/(H -E°) 



to the right-hand 



M = WMWi) [P/00 E%)] [9 - ((^-|(y + 22 - E^j)/(^j))] \iPj), 



(4k) 



which is the inhomogeneous bound-state analog of the Lippmann-Schwinger equation presented in Eq. (3e). 
The exact energy Ej of the exact physical bound state \if)j) is of course given by, 



Ej=({iPj\(V + H )\i>j)/(iPj\iPj)). 



(41) 
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The complete detailed prescription for the perturbational iteration of Eq. (4k) is then obviously, 



= I^X^f) - [P/(H [V ((^ n) \(V + H - E«)\^)/{^\^))] (4m) 

where, of course, = In addition, entirely as a byproduct of \ij}^ > ), we have, 

Ef +1) = ((Vi n) |(V + Ho)|^ n) >/^i n) |^ n) >) • (H 

As was pointed out earlier, the complete detailed prescription of Eqs. (4m) and (4n) for the perturbational 
iteration of a bound state approximation cannot be obtained within the confines of the standard Schro- 
dingcr bound-state perturbation approach [4] because that approach makes no attempt to work out the 
inhomogeneous Eq. (4j) or (4k) variant of the stationary-state Schrodinger equation for the exact bound state 
Of course having Eq. (4j) in hand is also absolutely critical for the development of nonpertubational 
"continued fraction" iteration of (pi\ipj), where (pi\ipj) must be nonzero for every member of the complete 
orthogonal basis set {{pi\}- The prescription for that nonperturbational "continued fraction" iteration clearly 
comes out to be, 

(Pi\4 n+1) ) = (Pi\^>j n) )/ 
{l + (M^y 1 (Pi\ [P/(Ho E]j\ [V ((^ n) \(V + H - E9)\^)/(^\^))} , (4 ° } 

where, of course, (pi\ipj°^) = (Pi\ipj)i an d the approximation Ej n+1 ^ to the energy eigenvalue is given by 
Eq. (4n). 

For a bound state approximation \ipj), finding a complete orthogonal basis set {{pi\} such that every 
{pi\Wj) is nonzero may not be simple because, aside from the ground state, bound states normally have nodes 
(where they of course vanish) in both coordinate and momentum representation. However if that bound 
state approximation is an eigenstate of an "unperturbed" harmonic oscillator Hamiltonian Ho, then 
one can find a related "harmonic-oscillator coherent false ground-state with false excitations" orthogonal 
basis {{ip(c)i\} such that every (^(c^lipj) is indeed nonvanishing. Here c is a complex number such that |c| 2 
is a transcendental positive real number, e.g., the natural- logarithm base e or the constant tt, and the "false 
ground state" is the minimum- uncertainty coherent state that is characterized by the complex number c, i.e., 
it is annihilated by one of the harmonic-oscillator annihilation operators with c subtracted from it. Successive 
mutually orthogonal states |^>(c)?) are then generated from this particular coherent "false ground state" 
characterized by c through the repeated action of the corresponding harmonic-oscillator creation operator 
with the complex conjugate c of c subtracted from it — these are the "false excitations" of the coherent 
"false ground state". It turns out that (^(c)?!^) can't vanish if |c| 2 is a positive transcendental number 
because it only vanishes when |c| 2 is a zero of an appropriate polynomial that has rational coefficients. We 
note that the nontranscendental (i.e., the "algebraic") real numbers, being a countable set, are of measure 
zero. Furthermore, choosing to use harmonic-oscillator states as bound state approximations is likely a 
viable proposition, e.g., for a bound state of a three-dimensional system, its approximation by a suitable 
three-dimensional harmonic oscillator state could probably be made to work out quite well. 

Time-dependent Dirac-picture successive approximation methods 

It is often the case that the Hamiltonian operator H for a time-dependent Schrodinger equation, 

ihd(\ip(t)))/dt = H\^(t)), (5a) 

can be written in the form H = Hq + V, where Hq is exactly diagonalizable and \ip(to)) is specified to equal 
one of the eigenstates | Vo) °f H - In that case it can be useful to reexpress the time-dependent Schrodinger 
Eq. (5a) in the Dirac picture [5], 

|V(*)> = e- iS °(*-*°)/ fi |V>B(*)} ; ( 5b ) 
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which yields, 

ihd(\i, D {t)))/dt = Vb(t)|Vz>(t)>, ( 5c ) 

where, 

V D (t) = f e +iHo(t-to)/hy e -iH (t-to)/h^ (5d) 

and V'-d(^o)} = NAo)- Through integration, Eq. (5c) can be reexpressed in an inhomogeneous form that 
incorporates the eigenstate |f/>o) of H as the inhomogeneous term, 

\il>o(t)) = |Vo) - / t * (5e) 

The inhomogeneous Eq. (5e) has the form of our generic Eq. (la), and therefore is subject to either standard 
perturbational iteration, 

IV4" +1) W> = lV>o> - {i/K)Jl V D (t')\^\t'))dt', (5f) 

where of course \i>£^(t)) = \ipo), or, if we can find a complete orthogonal basis set {{pi\} such that (pi\ipo) 
is never zero, we can carry out nonperturbational "continued- fraction" iteration, 



(5g) 



where of course (pi\ip^\t)) = (Pi\ipo)- 

A prime application of the Dirac-picture Eq. (5e) is to quantum field theories [6] , where it has invariably 
so far been pursued by using the standard perturbational iteration scheme of Eq. (5f), no doubt most 
famously by Dyson to systematically extract the quintessentially perturbational Feynman diagrams [6]. It 
would certainly be interesting if the nonperturbational "continued-fraction" iteration scheme of Eq. (5g) 
could conceivably be applied to quantum field theories. 

In quantum-field applications |^o) is a small number of unperturbed boson and fcrmion particles. Now 
unperturbed boson systems can be regarded as a massive collection of quantized simple harmonic oscillators, 
while unperturbed fermion systems can be regarded as an equally massive collection of quantum two-state 
entities. In the previous section we have pointed out that systems which are initially harmonic oscillators 
can be made compatible with continued-fraction iteration by using a "coherent false ground state" and 
related "false excitation" orthogonal basis (in quantum field theories perhaps better termed a "coherent 
false vacuum" and related "false Bose-particle" orthogonal basis), wherein quantized oscillator annihilation 
operators are translated by a complex number, and creation operators by that number's complex conjugate, 
which, of course can also be described as a simple real-valued c-number vector translation of the quantum 
coordinate-momentum operator phase space. The reason that such c-number operator translations simply 
generate additional orthogonal bases is, of course, that c-number translations have no effect on the crucial 
commutator algebra of those operators, i.e., they are canonical transformations, indeed unitary ones. 

Systems which are initially collections of two-state systems would seem in principle even easier to make 
compatible with "continued-fraction" iteration. In Pauli-matrix language, if the natural basis for the initial 
two-state system consists of the two eigenstates of |(1 + <r 2 ) (which has the two eigenvalues and 1, namely 
the two natural fermionic-state occupation numbers), then for this particular purpose of "continued fraction" 
iteration, we can instead adopt the two eigenstates of i(l + o- y ), for example, since each state of the natural 
basis indeed has a a nonzero inner product with each one of the states of this adopted basis. The example 
just given can be described as a spinor basis rotation of ninety degrees, and other rotation angles ought to 
achieve the desired result as well. Note that the natural unperturbed-particle basis has here once again been 
subjected to a unitary transformation. 

The fact that the continued-fraction iteration scheme forces on one a basis that has no zero-valued inner 
products with the members of the natural particle occupation-number basis of the quantum field system's 
unperturbed Hamiltonian operator H implies that this iteration scheme will very rapidly indeed force an 
unbounded number of unperturbed particles into active participation. That is in very sharp contrast with 
the perturbational iteration scheme: the intrinsic nature of Feynman diagrams immediately reveals that no 
more than a certain maximum number of unperturbed particles can ever participate through a given order of 
the perturbational iteration. (The only exception to this iron Feynman-diagram rule occurs when individual 
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diagrams fail by yielding infrared divergences - whose cure is widely agreed to involve coherent photon states, 
which have no upper bound on the possible number of photons which they describe). The continued-fraction 
iteration scheme, contrariwise, forces on one coherent states (i.e., "false vacua") and their related mutually 
orthogonal "false Bose-particle" brethren that all have an unbounded number of unperturbed particles. It 
as well forces on one the fact that all the fermionic basis states have a greater than zero probability to 
actually be occupied by unperturbed fermions . Feynman diagrams reveal "virtual particles", but continued- 
fraction iteration reveals immense "virtual clouds" which have no upper bound on the number of unperturbed 
particles involved. Landing in the thick of such "virtual clouds" within an iteration or two suggests much 
better convergence for the continued-fraction scheme than for the perturbational one. 

Feynman diagrams can diverge. In the continued-fraction scheme, everything that is calculated gets 
pushed into a denominator, which makes actual divergence quite difficult to visualize mathematically. What 
diverges in the continued-fraction scheme would seem to be unperturbed-particle participation rather than 
the iteration expressions themselves. 
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